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ABSTRACT 

We present detailed NLTE synthetic spectra of hydrodynamic SNe la 
models. We make no assumptions about the form of the spectrum at the 
inner boundary. We calculate both Chandrasekhar-mass deflagration models 
and sub-Chandrasekhar "helium detonators." Gamma-ray deposition is 
handled in a simple, accurate manner. We have parameterized the storage of 
energy that arises from the time dependent deposition of radioactive decay 
energy in a reasonable manner, that spans the expected range. We find 
that the Chandrasekhar-mass deflagration model W7 of Nomoto et al. shows 
good agreement with the observed spectra of SN 1992A and SN 1994D, 
particularly in the UV, where our models are expected to be most accurate. 
The sub-Chandrasekhar models do not reproduce the UV deficit observed in 
normal SNe la. They do bear some resemblance to sub-luminous SNe la, but 
the shape of the spectra (i.e. the colors) are opposite to that of the observed 
ones and the intermediate mass element lines such as Si II, and Ca II are 
extremely weak, which seems to be a generic difficulty of the models. Although 
the sub-Chandrasekhar models have a significant helium abundance (unlike 
Chandrasekhar-mass models) , helium lines are not prominent in the spectra near 
maximum light and thus do not act as a spectral signature for the progenitor. 

Subject headings: supernovae: general 
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1. Introduction 

Type la supernovae (SNe la) are among the brightest known objects in the universe. 
Since they form a nearly homogeneous class and simple selection criteria can make the 
observed dispersion quite small, they are natural cosmological probes (|Vaughan et al. 



1995|) . The observed homogeneity has led to a search for a homogeneous progenitor, that 
would satisfy the requirement of lacking hydrogen. This has led to the assumption that the 
SNe la progenitor involves the explosion of a Chandrasekhar-mass white dwarf. The current 



status of the search for the identification of the SNe la progenitor is reviewed in [Branch 



et al. (1995 ). Hydrodynamic explosion models have included deflagration models such as 



the "W7" model of [Nomoto, Thielemann, fc Yokoi (1984| ). While this model is somewhat 



hand-crafted to fit the observed spectra and suffers from an overproduction of neutron rich 
species, it remains the standard in the field. The "DD" (delayed detonation, [Khokhlov 



1991a|; Woosley 1991|) and "PDD" models (Pulsating Delayed Detonation, [Khokhlov 199TbD 



improve the predicted nucleosynthetic yield and gives qualitative agreement with the 
observed spectra and light curves ([Hoflich fc Khokhlov 199^ ). 



The observation of the super-luminous SN 1991T ([Jeffery et al. 19921) and the very 



sub-luminous SN 1991bg ([Filippenko et al. 1992| ) convincingly showed that the class of 
SNe la is not entirely homogeneous. In attempting to model SNe lb, Livne fc Glasner (1990 



|1991D examined a "helium-igniter" where a sub-Chandrasekhar or a Chandrasekhar-mass 
C/0 white dwarf with an accreted helium shell detonates near the center following the 
detonation of the helium shell. Motivated by observations of SN 1991T and 1991bg as well 
as other claims for evidence that SNe la form a sequence ( Phillips 1993|; Hamuy et al. 1995]) 



Woosley fc Weaver (199^) and p^ivne &: Arnett (1995| ) investigated the helium-igniter as a 



realistic model for SNe la. 

On the face of it helium-igniters have much to recommend them as a plausible 
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progenitor model of SNe la. By varying the initial white dwarf mass such models can 
produce a range of nickel mass M-^i 0.2 — 1.0 Mq, naturally leading to a sequence 
of supernovae. Sub-Chandrasekhar models do not suffer from the same neutronization 
problem that occurs in the C/0 deflagration models, and population synthesis studies may 
produce the progenitors in the requisite quantities (P?utukov, Yungelson, fc Iben 199^) , 
although not in old populations ( [Branch et al. 19951 ). 



On the negative side, the models produce an outer shell of the products of explosive 
helium burning: helium and ^^Ni; elements not typically associated with the outer layers of 
SNe la. In addition the light curves are extremely fast ( |Woosley fc Weaver 1994|; Hotiich~S 



Khokhlov 1996| ), so in particular, it is not clear that the observed photometric diversity can 



be reproduced by such models. 

Since one of the primary goals of synthetic spectral synthesis is the confrontation 
of theoretical models with observations, we present the results of synthetic spectrum 
calculations of the helium igniter models of [Woosley fc Weaver (199^ ) and |Livne fc Arnett 



|(1995|) . Since this is our first application of our program of synthetic spectral synthesis to 
hydrodynamical SNe la models, we also present the results of the spectral synthesis of the 
W7 model ([Nomoto, Thielemann, fc Yokoi 1984|) . 



2. Calculations 

The calculations are performed using the generalized stellar atmosphere program 
PHOENIX 7 . 1 ( IHauschildt 19924 b|, 1993|; Hauschildt, Baron, fc Allard 1996|) in generally 



the same way we have applied it previously to SNe la ( [Nugent et al. 19954b|) , although 



we have modified the code to allow the treatment of nebular boundary conditions, as well 
as stratified composition and a full gamma-ray deposition calculation. We have compared 
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the results of our 7-ray deposition with more detailed calculations and the agreement is 
excellent (Young & Kumagai 1995, private communication). The boundary conditions 
make no assumptions about the form of the flux at the inner boundary, but rather impose 
continuity requirements on the intensity (with correct Lorentz transformations). Thus all of 
the flux comes from the atmosphere itself, there is no "light bulb" at the center. PHOENIX 
accurately solves the fully relativistic radiation transport equation along with the non-LTE 
rate equations (for some ions) while ensuring radiative equilibrium (energy conservation). 
The following ions were treated in non-LTE in the calculations reported here (the number 
of levels follows in parenthesis): He I (11), He II (10), Na I (3), Ne I (26), Ca II (87), Mg II 
(18), C I (228), O I (36), Fe II (617), Co II (255), Ti II (204), S II (85) and Si II (94). 

The hydrodynamical models were evolved in time by assuming that the expansion is 
homologous, i.e. the velocity of any given mass point was held constant. The models were 
rezoned into 50 mass zones, with roughly a logarithmic spacing in Xstd, where Tstd is the 
total extinction optical depth in the continuum at 5000 A. Care was taken to resolve the 
density profiles. 

In theory, once we choose a time since explosion the model is completely determined. 
The density structure is specified by the homology transformation, the compositions are 
fixed (once decay of the radioactive species has been accounted for) and, since we use 
observed bolometric luminosity as an input parameter, the temperature structure is then 
completely determined by imposing the condition of radiative equilibrium. We parameterize 
the luminosity as: 

Lboi = vLf^ 

where L'^^ is the total instantaneous 7-ray luminosity deposited in the material and rj is 
a parameter that measures the net amount of energy stored over time by the material. 
Note that 77 differs from the parameter a defined by Arnett and co-workers ([Arnett 1982 
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Arnett, Branch, & Wheeler 1985) since a refers to the total instantaneous 7-ray himinosity, 
and instead corresponds to the parameter Q of [Hoflich fc Khokhlov (1996|) . [NB: While 



the definition of Q = « in [Hofiich fc Khokhlov (199"^ ), in previous papers in their series Q 



corresponds to r\\. Poflich fc Khokhlov (199^ found ri in the range 0.7 < r] < 1.8 for a wide 



variety of models that they examined, and we have varied i] from approximately 0.5 — 2.0. 
Actually 77 should be a function of radius, but an accurate calculation of rj will require a 
NLTE, multi-group radiation-hydrodynamical calculation. This procedure accounts for the 
time-dependent nature of the deposition of radioactive energy in an accurate manner. Given 
an input luminosity the temperature structure of the models is determined by demanding 
the modified radiative equilibrium condition: 

J ,^x{Bx-Jx)dX-S = 0, 

where S is the local instantaneous rate at which 7-ray energy is deposited. 



3. Results 
3.1. Model W7 

Figure displays our synthetic spectrum for the W7 model at 20 d past explosion for 
three choices of rj. This is several days after the time of bolometric maximum tboi = 14 d 
found by [Hdtiich fc Khokhlov (199"6|) and [Khokhlov, Miiller, fc H5tlich (1993|) who found 



= 1.3 at this time. The magnitudes and colors of these models are listed in Table |I|. The 
model strongly resembles observed SNe la spectra, showing the defining Si II A6355 line as 
well as the A5972 line, and lines from Ca II, S II, O I, and Fe II. Also, there is the strong 
UV deficit that is characteristic of SNe la spectra. The colors for the 77 = 1.0 20 d model 
are very similar (on average) to the colors found for normal SNe la. The = 2.0 colors are 
too blue and the rj = 0.8 are too red, suggesting that the value of r] = 1.3, found by Hofiich 
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and collaborators is reasonable. 

In addition to the NLTE lines that we treat directly, we must also include ~ 2 million 
additional lines in LTE. The shape of the spectrum is somewhat sensitive to the constant 
thermalization parameter e that we choose ( [Nugent et al. 1995a|; Baron et al. 1996| ), where 



e is defined by the source function for LTE metal lines, 

Si = il-e) I (t>uJvdp + eB,iT). 

Figure ^ compares the spectra for e = (10~^, 0.05, 0.1, 1.0). The UV (where most lines are 
treated in NLTE) is rather insensitive to the choice of e; however, in the optical, redward 
of 5000 A there is a strong dependence on e. Based on our previous work and the results 
of NLTE calculations, we choose e = 0.05 — 0.1 as our standard range. All the models 
discussed below have e = 0.1. 

Figure ^ shows the W7 model at 23 d (77 = 1.3) with an observed spectrum of SN 1992A 
[taken at 5 days after maximum light ( [Kirshner et al. 1993] )] and at 20 d (77 = 1.0) with 



an observed spectrum of SN 1994D [at maximum light in the optical and 3 days before 
maximum in the IR ( Meikle et al. 1996]) ]. The agreement is quite good across the entire 



range of observed wavelength for each supernova with all of the major (and most of the 
minor) features present in the synthetic spectra. While fine tuning could no doubt improve 
the fits, that is not our purpose in this paper. An interesting feature in Figure |^, is that 
both the observed spectrum for SN 1994D and the synthetic spectra of W7 show a "split" 
just blueward of the Ca II H&K feature. This is likely due to a blend of Ca II H&K and 
Si II A3858. [Kirshner et al. (1993| ) also noted that the two lines are of nearly equal strength. 



While the split is prominent in the observed spectrum of SN 1994D, it is clearly absent in 
the observed spectrum of SN 1992A. We will return to this issue in future work. 

Figure |^ shows the time evolution of the W7 model near maximum light. The 
bolometric magnitudes and colors of the 16 d, r/ = 1.1 and rj = 1.7 models (see Table 0) 



- 8 - 



should be compared with the results of [Hotlich fc Khokhlov (19961) who found Mboi = —19.56 
and B — V = 0.11 (for r] = 1.3). The bolometric magnitudes bracket the results of PofiicE 



|fc Khokhlov (1996|) , while both of the models are somewhat bluer than they found. This is 
likely due to differences in the treatment of radiation transport (NLTE, 82,000 wavelength 
points and full line profiles, versus LTE grey transport) and serves as an estimate of the 
theoretical uncertainty of such calculations. 



3.2. Sub-Chandrasekhar Models 



Figure || displays Model 2 of [Woosley fc Weaver (1994|) (WW2) 15 d after explosion for 



three choices of r]. This model is the explosion of an 0.7 Mq white dwarf that has accreted 
0.2 M0 of helium. Since these models are significantly less massive than W7, they peak 
earlier and hence it is sensible to examine them at earlier times. Figure ^ is similar to 
Figure |, but for Model 4 of [Livne fc Arnett (19951) (LA4), which is a 0.7 Mq white dwarf 



that has accreted 0.17 M0 of helium. The magnitude and color data for these models can 
be found in Table |^. The synthetic spectra have less line-blanketing and hence more flux in 
the UV, than does the W7 model. There is no strong evidence of either He I or He II lines. 
Although we do not use the Sobolev approximation at all in our calculations, we calculate 
the Sobolev optical depth of each NLTE line as a convenient diagnostic. The Sobolev 
optical depth of the He I A5876 lines is approximately 3 — 4 orders of magnitude weaker 
than that of the Si II A6355 line in both sets of models, thus these models are effective at 
"hiding helium." It is somewhat surprising that in a model with nearly 0.2 Mq helium on 
the outside that no evidence of helium should appear. This seems to be due to the very 
strong non-thermal ionization and the high UV flux which tends to keep the helium ionized, 
and/or highly excited, suppressing the strong He I lines. Optical He II lines, particularly 
Pa A4687.8 and Pp A3204.5, are also not prominent. At earlier times, with higher densities 
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recombination may populate the He I - II levels, but since the models will also be hotter at 
those times, it is not clear a priori that optical He I - II lines will ever be strong in these 
models. Understanding the exact suppression mechanism of helium lines will be the subject 
of future work. 

Figure |^ shows the synthetic spectra of the models of WW2 [i] = 1.1) and LA4 
{rj = 1.5) at 20 d, compared with the observed spectrum of SN 1994D. Clearly the synthetic 
spectra bear little resemblance to the observed SN la spectrum. The Si II A6355 line is 
quite weak and may not extend to high enough velocity. This is a generic problem with 
sub-Chandrasekhar models. The 15-d WW2 model of Figure ^ does however reproduce the 
boxy shape of the Ca IR triplet seen in SN 1994D, which may indicate that the calcium is 
confined to the correct velocity range in the model. 

Sub-Chandrasekhar models have been suggested as attractive models for low luminosity 
SNe la such as SN 1991bg. The spectra of this SN la [at maximum ( [t''ilippenko et aL| 
[I992D in the optical and lUE data in the UV ( [NASA Astrophysics Data I-'acilityTMll )] 
along with the spectra of the 20 d models of WW2 [rj = 1.1) and LA4 [rj = 1.5) can be 
seen in Fig. These models show some resemblance to the observed SN 1991bg spectrum, 
however, the shape of the spectra from these models is counter to that observed. This 
behavior is exemplified by the WW2 20 d model with rj = 1.1 in Table |^ where B — V is 
extremely red (0.74), but U — B is negative (-0.16). The fiat spectrum in the UV for these 
models (signified by a negative U — B) is observationally associated with bright SNe la 
such as SN 1991T and SN 1994D, rather than with dim supernovae such as SN 1991bg. In 
the spectrum of SN 1991bg, the trough near 4000 A is due to Ti II ( [^'ilippenko et al. 1992 



Nugent et al. 1995^1) . In the synthetic spectra this trough is not fiat enough to reproduce 
the observed spectrum of SN 1991bg, most likely because Ti II is confined to a small region 
in velocity in these models. Also, the calculated Si II and Ca II lines are very weak even 



- 10 - 



at this epoch. The strength of these hnes correlates inversely with the luminosity of the 
supernova ([Nugent et al. 1995b|) . Weak, dim supernova like SN 1991bg have strong Si II 
and Ca II lines, whereas powerful, bright SNe such as 1991T, have weak lines, thus a model 
that hopes to fit SN 1991bg should show prominent Si II and Ca II lines. 

Table |^ lists the relative concentration of the four most abundant species at three 
different velocities for each of the three models that we have examined (the models W7, 
WW2, and LA4 have i] = 1.0, 1.0, 1.5, respectively) at 20 days after explosion. The models 
have very similar compositions at low velocity, a mixture of highly ionized iron-peak 
elements and He II. The detonation models have somewhat more He II than does the 
deflagration W7. In the outer parts, the composition is different. Whereas W7 has singly 
and doubly ionized intermediate mass elements, the helium igniters have mixtures of neutral 
and singly ionized helium, and doubly ionized elements just below the iron-peak (calcium 
and titanium). At very high velocity {v > 20000 km s~^), the composition of W7 is mostly 
a mixture of C-0, while the WW2, and LA4 models are dominated by He I - II. As we 
have already noted the presence of intermediate mass elements at high velocity is required 
to fit the observed spectra. 

In order to further elucidate the differences between W7 and the helium-igniters. 
Figure || compares the temperature and electron density profiles of W7 {t] = 1.0) to WW2 
(r^ = 1.1) at 20 d. The more massive W7 model has a much steeper electron density profile, 
and it is cooler on the outside. Figure |TU| displays the Sobolev optical depth of the Co II 
-3 G° A2605.2 line as a function of r^td for the 3 models. While W7 and WW2 display 
similar Sobolev optical depths in this line at depth, LA4 only has a very small region where 
this line is optically thick, and, near the surface where the spectrum forms, both WW2 and 
LA4 are transparent in this line. This is typical for the iron-peak UV lines and it makes it 
very difficult for the helium-detonation models to display the proper line blanketing in the 



UV. It is not that the iron-peak elements are not present, but rather that they are not in 
the proper ionization/electronic states to create strong hne blanketing. 



4. Conclusions 

We have calculated very detailed NLTE synthetic spectra of hydrodynamical models 
for SNe la. We have used only symmetry considerations at the inner boundary and thus 
have not had to make any assumptions about the form of the flux there, the spectrum is 
calculated ab-initio. We have used a simple but accurate 7-ray transport algorithm and we 
have developed a reasonable parameterization of the time dependence of the 7-ray heating 
that can be compared with and calibrated to sophisticated radiation hydrodynamical 
calculations as they become available. 

The Chandrasekhar-mass deflagration model W7 shows good agreement with observed 
normal SNe la and it is likely that other Chandrasekhar-mass models such as DD or 
PDD ([Khokhlov, Miiller, fc Hofiich 19931) will also show reasonable agreement. While 



the sub-Chandrasekhar mass "helium igniter" models bear some resemblance to the 
sub-luminous SNe la typified by SN 1991bg the weakness of lines of the intermediate mass 
elements and the lack of the UV deficit will have to be addressed if these models are to 
remain as viable contenders for at least some SNe la. 
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Table 1. Magnitudes and Colors of W7 Models 



Time 


V 


MboI 


Mb 


Mv 


B-V 


U-B 


16 d 


0.6 


-18.48 


-18.64 


-18.60 


-0.04 


-0.35 


16 d 


0.9 


-18.85 


-18.95 


-18.78 


-0.18 


-0.64 


16 d 


1.1 


-19.10 


-19.11 


-18.92 


-0.19 


-0.79 


16 d 


1.7 


-19.61 


-19.44 


-19.26 


-0.17 


-0.91 


18 d 


0.9 


-18.70 


-18.93 


-18.83 


-0.10 


-0.37 


20 d 


0.8 


-18.51 


-18.71 


-18.81 


0.10 


-0.09 


20 d 


0.9 


-18.55 


-18.77 


-18.83 


0.07 


-0.12 


20 d 


1.0 


-18.69 


-18.93 


-18.93 


0.00 


-0.17 


20 d 


2.0 


-19.45 


-19.54 


-19.24 


-0.30 


-0.82 


23 d 


1.3 


-18.75 


-18.93 


-19.04 


0.10 


-0.14 


25 d 


1.5 


-18.81 


-19.11 


-19.08 


-0.04 


-0.20 



Note. — Mboi, Mb, and My are the bolometric, B and V absolute magnitudes of the 
models respectively. B — V and U — B are the associated colors. 
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Table 2. Magnitudes and Colors of the Sub-Chandrasekhar Models 



Model 


Time 


r) 


MboI 


Mb 


Mv 


B-V 


U-B 


WW2 


15 d 


0.7 


-18.08 


-18.04 


-18.68 


0.64 


-0.18 


WW2 


15 d 


1.0 


-18.49 


-18.72 


-18.82 


0.10 


-0.31 


WW2 


15 d 


1.3 


-18.78 


-19.13 


-18.83 


-0.30 


-0.46 


WW2 


15 d 


1.6 


-18.98 


-19.33 


-18.85 


-0.48 


-0.55 


WW2 


20 d 


1.1 


-18.22 


-18.21 


-18.95 


0.74 


-0.16 


WW2 


20 d 


1.5 


-18.51 


-18.78 


-18.97 


0.19 


-0.34 


WW2 


20 d 


1.9 


-18.80 


-19.17 


-18.82 


-0.36 


-0.65 


LA4 


15 d 


0.7 


-18.20 


-18.30 


-18.81 


0.52 


-0.17 


LA4 


15 d 


1.1 


-18.69 


-18.99 


-18.90 


-0.08 


-0.43 


LA4 


15 d 


1.7 


-19.17 


-19.41 


-18.91 


-0.50 


-0.76 


LA4 


20 d 


1.5 


-18.60 


-18.96 


-18.97 


0.01 


-0.39 


LA4 


20 d 


1.6 


-18.66 


-19.03 


-18.97 


-0.07 


-0.43 


LA4 


20 d 


2.0 


-18.89 


-19.29 


-18.91 


-0.37 


-0.60 



Note. — Mboi, Mb, and My are the bolometric, B and V absolute magnitudes of the 
models respectively. B — V and U — B are the associated colors. 
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Table 3. 4 Most Abundant Species for Each Model at 20 Days After Explosion 



Model 


Velocity 




Relative Concentration 






(km s~^) 




(ppm) 




W7 


20,000 


C I (2 X 10^) 


II (2 X 10^) 


Ne II (2 X 10^) 


I (1 X 10^) 


LA4 


20,000 


He II (4 X 10^) 


He I (9 X 10^) 


Ti III (6 X 10^) 


Ca III (4 X 10^) 


WW2 


20,000 


He II (4 X 10^) 


He I (2 X 10^) 


Ti III (7 X 10=^) 


Ca III (6 X 10=^) 


W7 


10,000 


Si III (1 X 10^) 


Co III (7 X 10^) 


S III (6 X 10^) 


Ca III (3 X 10^) 


LA4 


10,000 


Si III (1 X 10^) 


Co III (1 X 10^) 


S III (5 X 10^) 


Ni III (1 X 10^) 


WW2 


10,000 


Si III (2 X 10^) 


S III (1 X 10^) 


Ar II (1 X 10^) 


Ca III (1 X 10^) 


W7 


5,000 


Co IV (2 X 10^) 


Ni IV (7 X 10^) 


Fe IV (2 X 10^) 


He II (2 X lO^) 


LA4 


5,000 


Co IV (2 X 10^) 


He II (7 X 10^) 


Ni IV (2 X 10^) 


Fe IV (2 X 10^) 


WW2 


5,000 


Co IV (2 X 10^) 


He II (3 X 10^) 


Ni IV (2 X 10^) 


Fe IV (2 X 10^) 



Note. — The relative concentrations by number of the four most abundant species (in 
parts per million) for a particular zone (labeled by its velocity). The models W7, WW2, and 
LA4 have r] = 1.0, 1.1, 1.4, respectively. 
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Fig. 1. — The synthetic spectrum for the W7 model at 20 d past explosion for three choices 
of r). 



-18- 




Fig. 2. — The synthetic spectrum for the W7 model at 23 d past explosion (77 = 1.3) for 
four choices of e. The model with e = 0.10 is in radiative equilibrium and the temperature 
structure has been held fixed at that structure for all the models displayed. 
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Fig. 3. — The W7 model at 23 d (r^ = 1.3) compared with an observed spectrum of SN 1992A 
[5 days after maximum hght ( [Kirshner et al. 1993| )] and at 20 d (77 = 1.0) compared with an 
observed spectrum of SN 1994D [at maximum hght ( [Meikie et al. 199^ ]. 
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Fig. 4. — The time evolution of the W7 model near maximum light (ry = 1.1, 0.9, 1.0, 1.3, 1.5 
at t = 16, 18, 20, 23, 25 d, respectively). 
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Fig. 5.- 



Model WW2 15 d after explosion for three choices of rj. 
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Fig. 6.- 



Model LA4 15 d after explosion for three choices of rj. 
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Fig. 7. — The synthetic spectra of WW2 [rj = 1.1) and LA4 [r] = 1.5) at 20 d compared with 
the observed spectrum of SN 1992A [5 days after maximum hght ( [Kirshner et al. 1993| )]. 
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Fig. 8. — The synthetic spectra of WW2 (77 = 1.1) and LA4 (77 = 1.5) at 20 d compared 
with the observed spectrum of SN 1991bg [at maximum hght ([Fihppenko et ah 1992|; NASA] 
Astrophysics Data Facihty 1991| )]. 
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Fig. 9. — The electron density and temperature profiles for WW2 and W7 at 20 d, with 
r] = 1.1, 1.0 respectively. The "ghtches" in the density and temperature occur at the edge of 
the burning fronts and represent real changes in model structure. 
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Fig. 10.— The Sobolev optical depth of the Co II G° A2605.2 line as a function of 

Tstd for the 3 models listed in Table |[ The solid line denotes the W7 model, the dashed line 
WW2, and the dot-dashed line LA4. Note that rsoboicv is a purely local quantity and hence 
is not necessarily monotonic. 



